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Abstract 


We compute the transition amplitndes between charged particles 
of mass M and m accelerated by a constant electric field and inter¬ 
acting by the exchange of qnanta of a third field. We work in second 
qnantization in order to take into account both recoil effects induced 
by transitions and the vacunm instability of the charged fields. In 
spite of both effects, when the exchanged particle is neutral, the equi¬ 
librium ratio of the populations is simply exp(7r(M^ — vn?‘)/eE). Thus, 
in the limit (M — m)/M 0, one recovers Unruh’s result character¬ 
ized by the temperature a/27r where a is the acceleration. When the 
exchanged particle is charged, its vacnum instability prevents a sim¬ 
ple description of the eqnilibrinm state. However, in the limit wherein 
the charge of the exchanged particle tends to zero, the eqnilibrinm 
distribution is once more Boltzmanian, bnt characterized not only by 
a temperature but also by the electric potential felt by the exchanged 
particle. This work therefore confirms that thermodynamics in the 
presence of horizons does not rely on a semi-classical treatment. The 
relationship with thermodymanics of charged black holes is stressed. 
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1 Introduction 


Shortly after Hawking’s seminal discovery of black hole radiation [|I|, Unrnh|^ 
showed that it possesses a flat space analogue, namely that a uniformly ac¬ 
celerated detector perceives Minkowski vacuum to be thermaly populated at 
temperature Tu = a/2Ti. In Unruh’s original work, only the detector’s in¬ 
ternal states where treated quantum mechanically. Its position was treated 
classically and thus insensitive to the transitions. This is clearly an approx¬ 
imation since it violates momentum conservation. In order to analyze the 
validity of this approximation and to re-establish momentum conservation, 
one must treat the detector position quantum mechanically [^, |^. This 

enlargement of the quantum dynamics to a variable formally external to the 
dynamics provides new insights about the Unruh process. Moreover, it may 
also be used as a guide to other problems dealing with particle creation in 
the presence of horizons since, in those cases, similar approximation schemes 
are used and should be abandoned in other to address the question of the 
quantum back reaction. The new insights concerning the Unruh effect which 
have been obtained are: 

1. The detector can be described by a delocalised wave function, where¬ 
upon the classical geometric notion of a horizon no longer exists. (Of 
course, one may approximatively recover the concept of horizon upon 
building well localized wave packets.) Nevertheless thermal rates for 
transitions of the detector still obtain thereby conhrming that thermo¬ 
dynamical relations still govern the physics when one goes beyond the 
semiclassical treatment. 

2. Each time the detector makes a transition, it recoils both in momentum 
and in energy in such a way that the total instantaneous Minkowski 
momentum and energy is conserved in the sense explained in |^, pages 
245-246. 

3. These recoils give rise to a decoherence of the detector-radiation sys¬ 
tem. This in turn implies that the detector emits a steady flux of 
radiation, contrary to the situation where the detector’s position is 
treated classically ||^, |^, ||. 

4. The Unruh effect and the Schwinger process[^, of pair creation in an 
electric held are “in equilibrium” . Moreover the area of the acceler¬ 
ation horizon plays the role of an entropy in delivering the equilibrium 
population ratios. 
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To understand the dynamical origin of these equilibrium notions is one of our 
main concerns. Thus, we hrst present the content of point 4 in some detail. 
In [^, the accelerated detector is modeled by a “two level ion” propagating 
in constant electric held E. It has charge Q and its two levels have rest mass 
M and m. It therefore uniformly accelerates with acceleration QE/M or 
QE/m according to its mass. The ion can make transitions between its two 
levels by emitting or absorbing a massless quantum. Thus it behaves like an 
accelerated particle detector with mass gap AM = M — m. Moreover, when 
describing the ion’s states by held operators, i.e. by working in second quan¬ 
tization, the electric held leads to vacuum instability through pair creation 
of ions anti-ions. The mean numbers of created pairs of mass m and M are : 


A priori independently of these creation processes, an ion will make transi¬ 
tions from its excited to its ground state at a rate or from its ground 

to its excited state at a rate Rm^M- As suggested in [O and conhrmed in 


, these transitions are intimately related to the creation processes since 


the ratio of their rates is given by 


Rm^M _ Nm _ ^-TT{M^-m?)/QE 


R 


M^m 




( 2 ) 


Thus we can say that the two processes are in equilibrium since they deter¬ 
mine the same distribution. We emphasize that eq. is exact in the sense 
that it takes into account all ehects due to the hnite mass of the detector, 
i.e. recoil ehects, and the hnite probability to create pairs of detectors. 

Upon taking the limit M,m —>■ oo, with a = 2QE/{M -|- m) and AM 
constant, recoil ehects and pair creation amplitudes vanish. Therefore, one 
expects to recover Unruh’s result which gives the equilibrium probabilities of 
an accelerated detector of “inhnite” mass and of given acceleration a. Indeed, 
upon taking the above limit, eq. (|) becomes 


Rn 


R 


_ ^-TTAM{M+m)/QE _ ^-2vAM/a 


M^m 


( 3 ) 


The interest of eq. (H) is further enhanced when one recalls that the 
probability for pair creation due to the Schwinger process can be expressed 
in terms of the change of the area of the acceleration horizon [|I^ p!4| 


7Y — 


( 4 ) 


where AAnifn, Q) = — Auim, q) is the area of the acceleration horizon 

in hat space minus the area if a pair of particles of mass m and charge e is 
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emitted^ Thus eq. (^) can be rewritten as 

^ ^-iAAH(M,Q)-AAHim,Q)} ^ ^-AAh{AM) 

RlVl^m 

where in the last equality we have written AAh{AM) = Anim, Q)—Ah{M, Q) 
as the difference of the horizon area between the initial and final states. The 
area of the acceleration horizon therefore behaves like an entropy in delivering 
the equilibrium population ratios. 

The aim of the present paper is to extend these notions to the situation 
wherein the field <h with which the detector interacts is charged and therefore 
feels also the constant T^-field. Our main results are: 


1. The transition rates and the equilibrium distribution are no longer 
described by Boltzman ratios. Hence the connection with thermody¬ 
namics no longer obtains. Similarly the equilibrium relation, eq. (^) 
between Schwinger and radiative processes is broken. 


2. However in the case of weakly charged exchanged particles, i.e. for 
a/Q where a is the charge of the field <h, an extended thermody¬ 
namical relation is obtained. Indeed, in that case, eq. becomes 


^ )+0(aVQh ^ (5) 

In a forthcoming article we shall confirm these results by adopting the more 
kinematical point of view of quantizing d) in Rindler coordinates ||15||. 

The thermodynamical relation eq. (j^) shows that, in addition to the 
Unruh temperature Tjj = a/2Ti, there is now an electric potential (= E/2a) 
which modifies the equilibrium. This is strictly analogous to the ratio of the 
rates for charged particles to be emitted or absorbed by a charged black hole. 
In Hawking’s derivation, this ratio is expressed in terms of the Bogoljubov 
coefficients 7 ^^, (3^ characterizing the mixture of in and out modes of the <I> 
field. However, to make clearer the contact with eq. (j^), we can express 
Hawking’s result in terms of the rate IZm^m-uj to jump from a black hole of 
mass M and charge Q to a black hole characterized by M — a; and Q — a 
and the rate of the inverse process Indeed, Hawking’s result can 

be expressed as 


1 2 _ ^l3H{uj—aip) _ 3Zm—u)^M 

Juj TZm^M-lo 


(7) 


® More specifically it is the change of the area of the horizon in the euclidean continu¬ 
ation of space time, see the cited references for additional details 
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where is the inverse Hawking temperature, lo is the energy of the quan¬ 
tum measured at spatial inhnity and 0 is the difference of electric potential 
between the horizon and inhnity. In eq. (^, the equivalent of (j) is, E/2a, 
the difference of electric potential between the horizon and the accelerated 
trajectory where the charged quantum is emitted (absorbed). Indeed, in 
Rindler coordinates r, p, the (static) potential is 


A^{p) = Eap^/2 


( 8 ) 


and the particle is localized at p = 1/a. 

We wish to emphasize that in both cases the thermodynamical canonical 
relations, eqs. ( 0 . 0 ). can be a posteriori related to changes in horizon’s 
area. Indeed, in both cases, the logarithm of these ratios is equal to the a 
quarter of the horizon area change when linearized in the energy and charge 
differences, i.e. in the energy and charge carried by the quantum of the $ 
held. In the black hole case, this is a reexpression of the hrst law of black hole 
thermodynamics. In the present case, since eq. (^) deals with the change in 
the accelerated horizon. 


AAffiAM^a) _ ^E^ 

. —linearized _ 

4 a 2a 


(9) 


should be considered as the (linearized) hrst law of accelerated-horizon ther¬ 
modynamics. 

However, in contradistinction to the treatment used to derive eq. (0) , we 
shall not make use of any background held approximation to obtain eq. (^. 
Instead, we shall approximate the exact result by taking variations limited to 
hrst order in a and AM, see Section 4. Therefore, we can analyze the hnite 
diherence and not only the hrst order changes delivering the above canonical 
concepts through diherentiation. Upon neglecting the vacuum instability 
of the exchanged $ held, we obtain that the logarithm of the ratio of the 
transition rates is given by the following hnite diherence 


AH 


H 


= vrl 


~ qE^ 


( 10 ) 


where q = Q — a is the “residual” charge of the ground state of mass m. 
Since eq. ([T 0 |) agrees with the Schwinger process, see eq. (D, by virtue of 
eq. (^, this hnite diherence is equal to the hnite change in area, as in eq. 
d^). This conhrms that a quarter of the area not only delivers canonical 
distributions and thermodynamics but truly determines quantum processes 
as in statistical mechanics. 
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To conclude let us recall once more the origin of the agreement of eq. 
d^) with Unruh’s result (or the agreement of eq. (§) with eq. (|^)). Eq. 
(|^) is derived in an enlarged quantum setting wherein the trajectory of the 
heavy ion is quantum mechanically treated whereas Unruh’s derivation is 
based on a background held approximation since the trajectory is classically 
determined. The agreement of the ratio of the transition rates evaluated 
from both treatments arises when the following procedure is applied to tran¬ 
sition amplitudes evaluated in the more quantum framework. Upon working 
with WKB waves and performing hrst order expansions in AM/M and in 
the momentum transfer, these transition amplitudes coincide with the cor¬ 
responding amplitudes evaluated at the background held approximation, see 
for the details. What guarantees this agreement is that hrst order expan¬ 
sions of the WKB phases are controlled by Hamilton-Jacobi equations. The 
same relation will therefore hold when one considers gravity. (This has been 
explicitly verihed in mini-superspace in [|^].) To hrst order in the matter en¬ 
ergy change, transition amplitudes computed in quantum gravity with WKB 
waves are equal to the corresponding amplitudes evaluated from quantum 
held theory in a given classical geometry. Only second order changes, i.e. 
the non-linear response of gravity, involve the Planck mass. 

The present article is organized as follows. Section 1 is this introduc¬ 
tion. Section 2 is devoted to recalling the quantization of a charged held 
in an electric held and the Schwinger process. In Section 3 we introduce 
our detector model, and consider the case when the detector interacts with 
a massive, but neutral, held. The techniques developed in this section are 
then used in Section 4 to analyze the more complicated case of interactions 
with a charged held. Appendix 1 is devoted to the analytical evaluation of 
the transition amplitudes and Appendix 2 deals with the limit a i—0 of the 
transition amplitudes. 


2 Charged Particles in an Electric Field 


The aim of this section is to review the quantization of a massive charged 
scalar held in an external electric held E. For the reader interested in a more 
complete treatment we refer to and references therein. Classically, the 
equation of motion of a relativistic charged particle of mass M and charge 
Q in an electric held is: 


M 


(Px^ 

dP 


-QE 


dXy 

dr 


( 11 ) 


5 






Its trajectory is a hyperbola with parametric equations given by: 


X — Xq = — cosh ^ 4(5 — So) 


( 12 ) 

(13) 


Here A = QE/M is the classical acceleration of the particle The time 
coordinate of the turning point of the hyperbola is to = —k/QE where 
k = {Mx + Ax) is the conserved momentum canonically conjugate to the 
variable x in the gauge At = Q and Ax = —Et. 

The corresponding Klein-Gordon equation for the held is: 


□ + M^ '4)M{t,x)= {-id, 


+ QEtf + dl + M^l x) = 0 (14) 


The general solution of this equation can be written as a superposition of 
modes: 

^ikx 

rR — X.M,k 
y I'X 


(t) , 


(15) 


V’M,fc(t,x) 


and their complex conjugates. The functions XM,k{t) obey the equation: 


dt + + (fc + QEt)^ XM,ki^) ~ 0 


(16) 


Because of the t dependence of this potential there will be some backscat- 
tering. To interpret this in the context of the Klein-Gordon equation, we 
note that the current operator is Jq = —iQ d f Thus the backscattering 
corresponds to mixing of positive and negative charge. In a second quan¬ 
tized context this corresponds to pair creation 0.111. To describe them we 
introduce a set of solutions with only positive or negative charge for t —»• —cx) 
(in-modes) 


rzit^x) 

^M?-k{t^x) 


^ikx g-ftM 
v/^(2QE)V4^---V2 
pikx 

p—ikx 




( 17 ) 

(18) 


and a set of solutions with only positive or negative charge for t —>■ -|-cx) 
(out-modes) 
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Akx 


La out 


[ X] 

(19) 

PuJ-ki-t^x) 


p—ikx 


^v* \ x\ 

(20) 


where eM = M^/ {2QE) and where the superscripts p and a refer respectively 
to particle and anti-particle wave functions. We have introduced a synthetic 
notation for the parabolic cylinder functions and their argument 


= 


(2QE)V4 

k 


D 








( 21 ) 


We also note that the parabolic cylinder function have the following integral 
representation |T7[| : 




g 2 g 27r/8 


,.i \2 

-e 2 


« POO _ _ ,,2 1 

X^QE / ^-iX^2QEv+i — ^-ieM -2 (22) 

Jo 


{2QEy/^T{l-teM) 

where the integration parameter v is classically related to t and its conjugate 


momentum pt by|18 


V = jQE/2 


Pt / kn . 



,As 


Since the parabolic cylinder functions have the property: 

= aMV.„i-\]+PMv:j-x] 

the in- and out- modes are related by the linear transformation: 

where 


(23) 


(24) 


CtM = 


\/^e * 


* 4 g 2 


r(| -|- Km) 


and 


Pm = ie 




(25) 

(26) 


These Bogoljubov coefficients are k independent, but mass and charge de¬ 
pendent. The second quantized held Tm should be decomposed in both the 
in or out basis 


TM(t, x) = dk 


>{ 




(t,^) (k) + (PmT ipx))* (k) 


I 




, (27) 
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to define the in and ont operators. From eq. (p5D we obtain 


aT = aMaT^-PMhTl\ 


( 28 ) 


The Heisenberg state | 0, in) contains no particles at early times, i.e. it is 
annihilated by the in destrnction operators. At late times it contains pairs 
of particles, as expressed by the relation 

, , Pm r, Olitf J Olitj T7 

|0,m)M = \0,out}M (29) 

where 

= m( 0 , out I 0, in)M- (30) 

is a normalization factor. The mean nnmber of created particles is 


3 Particle interactions: emission of a neutral 
particle 


In this section we consider a nniformly accelerated detector interacting with 
a nentral scalar held of mass fx. We shall show that eq. (0), i.e. the 
relation between Schwinger and radiative processes which was obtained in 
for a massless held, still holds when the exchanged qnanta are massive. 
Moreover, we shall see that all amplitudes linear in the coupling constant can 
be expressed in terms of a single amplitude describing the creation of a pair of 
charged quanta. The techniques developed in this section will be generalized 
in the next section where the held is both massive and charged. 

The detector is described by two charged scalar helds and with 
mass M and m and the same charge Q, propagating in the electric held E. 
As in O, H, O , the interaction between the helds is supposed to be given by: 


j^int _ 


= 9 


dx 


+ <h,. 


(32) 


A hrst amplitude of interest is the amplitude A of transition from an in 
M-particle of momentum k into an out m-particle of momentum k' and 
a /r-particle of momentum k” , which we schematically write as M(/c) —> 

m{k') + jjiik"). It corresponds to spontaneous deexcitation of the detector 





(since M > m). In the interaction representation, to first order in g, it is 
given by 


A{k\k\k") 


—ig [ dtdx{0,out \ a^{k'') a°^^{k') I 0,m) 

(33) 


The factorization of the vacuum states as tensor product of three vacua: 


leads to the expression (see [|T^ ^): 

A{k\k\k'') = ^ ^ {-i)g f dtdx 


(34) 


(35) 


where A/m, A/’m are the overlaps (see eqs. (|^J30|) ). 

At this point it is intersting to note that the quantity A{k\k\ k") describes 
several different processes in addition to eq. (^). Indeed it is not difficult 
to verify using eqs. (|T7|-^0|) that the four transitions 

M{k) m(fc') + fi{k'') 

M{-k) m{-k') + /x(-fc") 

m{k') + /i(fc") ^ M{k) 

m{-k') + /r(-A)") ^ M{-k) , (36) 


where M and ffi denote antiparticles all have the same amplitude A{k\k', k"). 
These different amplitudes are related by combinations of substituting par¬ 
ticles for antiparticles, changing the sign of the momentum k ^ —k and 
permutting incoming and outgoing quanta. Their origin will be further ex¬ 
plained in the next section. Similar properties will obtain for all the am¬ 
plitudes we shall introduce. Thus the notation we use for A{k\k',k"), and 
which is also used for all the other amplitudes in this paper, is that the | 
separates incoming from outgoing quanta, and that is the momentum of 
the Tm quanta, ±k' of the quanta, and ±/c" of the <I)„ quanta. 

Similarly, the amplitude B of spontaneous excitation of an m particle into 
a M particle accompanied by the emission of a /i particle is given by 


B{k'\k,-k") = -ig J dtdx{0,out I aA-k") a°A\k) a“t(fc') | 0,*n) 
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Due to the uniform acceleration, this spontaneous excitation amplitude 
is non vanishing. As in the Unruh treatment, the ratio of the rates of spon¬ 
taneous excitation to spontaneous deexcitation is simply given by 

Rm^M _ |gp 

Riviera 


(38) 


since the norm of the amplitudes are independent of the momenta k and k'. 

Thus, when the detector reaches equilibrium, the ratio of the probabilities 
to hnd the detector in its excited or ground state are 


P, 


M equil 


R. 


■m^M 


P, 


m equil 


R 


M^m 


1 ^1^ 


(39) 


Our task is to calculate this ratio and to conhrm the relation with the 
Schwinger process, eq. (^Tj). To this end we introduce a third amplitude 
V corresponding to the creation from vacuum of an out-M-anti-particle, an 
out-m-particle and a /i-particle. To hrst order in g, it is given by 


V(| — k,k',k”) = —ig j dtdx{0,out \ a^{k 

1 1 


C*(fc') bZ\-k) I 0,tn) 


Rfm 


-ig) 


dtdx 


(40) 


Using the Bogoljubov transformation for the M-particle wave function, eq 
(1^) , we can reexpress V in terms of the A amplitude 


V{\-k,k'A") 


-w 


dt dx 


au-rM'"-: + Ml"! 


Mr 

[3MA{k\k', k") aM^ik, k', k 


p out 



(41) 


Similarly, using the Bogoljubov transformation for the m-particle we can 
reexpress V in terms of the B amplitude 

V{\k,P,kl = PmB{k'\k,-k") + arr,P{k,k',k") . (42) 

Using the identities (which constitute the main mathematical result of this 
section, and are proved hereafter) : 

I = T = 0 (43) 


one obtains 


(3mB{k'\P -k") = V{\-P k\ k") = l3MA{k\k\ k") , (44) 
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figure 1 


Figure 1: A graphic representation of the relation between amplitudes B, V, 
A of section 3. Particles are denoted by M{k), m{k'), fi{k'') according to 
their mass and momentum, and antiparticles by M, rh (their is no antiparti¬ 
cle for the held which is real). The momenta in each diagram are always 
conserved k = k' + k". The coefficients Pm, Pm which weight the different 
amplitudes are the amplitudes to produce a pair of M, M or m, m particles 
from vacuum. The picture is oriented both in time and space: particles 
(antiparticles) are incoming from, and outgoing to, the right (left) accord¬ 
ing to their acceleration in the electric held. The chargless particles are 
represented as vertical lines because they are not accelerated. 

relations that we have depicted schematically in hgure 1 . Inserting eq. (|4^) 
into eq.(^) yields the link between radiative and Schwinger processes dis¬ 
cussed in the introduction eq. (|) 

Pm eqmi _Nm _ 

Pm equil 

It remains to establish eq.(^). To this end, we write 

I(k,k',k") = Cst fdxdte‘’“e-‘’‘'‘e-‘’‘"^V,„l-\]V;j\']e‘“"‘ 

= 2iTS{k -Id- k") I dW,J-\]V:j\’]e‘"“‘ ( 46 ) 

where we have introduced X = {t + k/QE), X' = {t + k'/QE). To per¬ 
form the t integral we replace the Whittaker’s functions TX by their integral 
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representations eq. 


to obtain 


X(k, k', k") = S(k -k'- k") / e“"‘ 


(47) 


dudv e 


+iuy/2QE{t+^E)+iv?+i{t+^E)y/2QEv-iv'^/2 iem-l/2 
LL 


The quadratic phases in t cancel, and the t integral yields + n + )) • 

The argument of the delta function never vanishes on the domain of inte¬ 
gration of u and v since ut" + k" ^ 0 for /i 7 ^ 0. Hence X = 0. A similar 
reasoning shows that X' = 0. 

The interested reader will now hnd the calculation of the amplitudes 
themselves. By virtue of eq. (^), we only need to calculate the amplitude 
V. We start from eq.(@ in which we reexpress + 

The first term does not contribute since it is equal to 
X(fc, k', fc")/|ap. The second term gives: 


V(|-fc,fc',fc") = Cst dxdte'^'^^e-^'^'^e-^^''^V,^[X\V*JXy‘^''* 


= Cst' 5 {k - k'- k") j dt e Ef 

[ dudv g-™V^QE(t+-^E)+iu^/2^-'‘M-k/2 ^i(^t+^E),y2QEv-iv'^/2^iem-l/2 


>0 


pOO 

= Cst''S{k — k'— k'') / dudv S{—uV-\- 

Jo 


u 


-)) 


/O VWE 

^i{^E)y/^v-zvC2^ie,n-l/2^-"M-k/2^-^uy/^i^E)+iuC2 

roo . uj-k" f ^ 

= Cst'" 5 {k-k'-k") dv e'" 

Jo 

The last integral gives an integral representation of a Whittaker’s function 
(See [^, formula 3.383.4). Reinstating the value of Cst"' yields 


(48) 


Vi\-k,k',k")=g 


1 


aM-XlMOim-Mn 


-6{k -k' - k") 


StT 377 

e~ “3" ^ e~"3" 

2QEV^ 


(49) 


\LO+k"J 


— f (em+eM)!i(^m —Cm) 




2QE 


It is interesting to note that the delta function d(—u -|- n -|- ^^ 7 ^^) has the 
interpretation of ensuring local energy conservation. Indeed, replacing u and 
V by their classical relation to t and pt eq. (|^) , yields simply —= 0 
where pt is the energy of the M-particle, p( of the m-particle, and uj of the p 
particle. 

We conclude this section by calculating the p —> 0 limit of the above 
amplitude. This will allow us to make contact with the expressions of 
that were obtained using a different method. The evaluation of this limit 


12 




















needs some precautions. Indeed for k" < 0, the support of the delta function 
given by the integral over t in eq.(^) rejoins the boundary of the domain of 
integration. In order to avoid ambiguities, we have to substitute in eg. (p}]) , 
according to the sign of k", the decomposition of or in 

terms of in and out helds and take into account the inhnitesimal imaginary 
part that the squared masses of the helds share. For small value of /i one 
hnds: 




1 g-f (em-EM) g3j7I-/4 

’ 2QEV^ VmE 

Oh" 1 

QE r(^ + ieM)r(l + *em ~ 




Oh" 

QE^ 


r(^ + m )r(l-*e 

m + if-n) 



Note the surprising fact that in this limit, the amplitudes of decay of the same 
process but with the opposite momenta differ. However, this is peculiar to 
two dimensions. Indeed in higher dimension, the mass /r contains the squared 
transversal momentum and so even in the massless limit, vanishes only on a 
domain of zero measure in phase space. 


4 Particle interactions: emission of a charged 
particle 

In this section we shall suppose that the held is also charged. Thus there 
are three helds '^m, with masses M, m, fi and charges Q, q, a 
respectively. These three helds interact through the hamiltonian 


= g J dx 


(51) 

Charge conservation requires 




Q = q + a 


(52) 

The amplitude A of transition from an in M-particle into an out m-particle 

and a out /i-particle is given 

at hrst order in perturbation theory by: 


A{k\k',k") = {0,out\ 

Cik") a-\k'){- 

■i) J dtW’^AZ\k) 1 0,m) 


1 

m OinJsI ^ 

1 jj. J ipout 1 pin* 1 pin* 

J dtdx Cfc' r^k" 


= C6{k- 

k’ -k’’)A{k\k’,k”] 

1 

(53) 
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where 


1 


C = 


-w 


y/‘2,T^ CKm-^ 

and we have introduced the reduced amplitude 

r+oo 


/ -l-OO 

-CO 


(54) 


(55) 


and used the abbreviated notation for the parabolic cylinder functions in¬ 
troduced in eq. (pi]). Note also the appearance of the factor 1/a^Af^ which 
arises due to the vacuum instability of the fi held. 

As in the previous section, the processes 


M(fc) —> m{k') + , 

M{-k) m(-fc') + fl{-k'>) 

m{k') + fl{k") M{k) , 

m{-k') + fi{-k”) M{-k) , (56) 


all have the same amplitude A{k\k', k"). Note that these differ slightly from 
eq. (^) because one has to distinguish between /i and /2 since the held 
is charged. 

Similarly the amplitude B for a anti-m-particle to spontaneously excite 
into a anti-M-particle and a /i-particle is 

Bik'\k, -k") = -^ I dtdx{0,out I a°;\k") hZ\-k) \ 0,zn) 

= C6{k-k' -k”)B{k'\k,-k”) (57) 

where 

/ + 00 

( 58 ) 

-CO 

Another amplitude that we need to introduce is the amplitude C for a anti- 
/x-particle to spontaneously excite into a anti-M-particle and a m-particle: 

C{k"\ - k',k) = -ij dtdx{0,out I b°^\-k') aZ\k) I 0,m) 

= CS{k-k' -k”)C{k”\k,-k') (59) 


where 


C{k''\k,-k') 





(60) 
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and the amplitude V for spontaneous creation from vacuum of an anti-M- 
particle, a m-particle and a /i-particle 


V(| - k, k\k") 


-i I dtdx{0,out I | 0,tn) 

CS{k-k’ -k’')V{\-k,k’,k’') (61) 


where 


/ + 00 

( 62 ) 

-oo 

The rules governing the product of T> functions in the reduced amplitudes 
A, B, C, V can be summarized as: 

• The sum of the momenta is conserved. We always assume that k = 
k' + k”. 


• Equations (|^ and tell us that a function A] is associated to 

an incoming particle (of mass M, charge Q > 0, and momentum k) or 
an incoming antiparticle (of mass M, charge —Q, and momentum —k) 

• Equations (0 and 1^) tell us that a function V*^ [A] is associated to 
an outgoing particle (of mass M, charge Q > 0, and momentum k) or 
an outgoing antiparticle (of mass M, charge —Q, and momentum —k) 

• The change of variable t —>■ —t in the integrals accompagnied with a 
change in sign of momentum k —k replaces incoming quanta by 
outgoing quanta and vice versa, without changing the value of the 
integral. 

These rules give directly the relation eq. (|56D and similar ones for B, C and 
V. Therefore the 64 possible amplitudes describing Erst order interactions 
between and quanta can all equal to either A, B, C, VQ. 

There is a further identity which will play a crucial role, namely 

/ + 00 

[A"] = 0 . (63) 

-OO 

It generalizes eq. (|43|) and it is proven also by using the integral representa¬ 
tion of the T) functions given in eq. (^^j). As in eq. (E3). one verihes that X 

^Note that these amplitudes are all integrals of products of three V* functions. Math¬ 
ematically one could also consider intergals of products of three V* and V functions. The 
Bogoljubov transformation ( p 5 | ) ensures that they can be reduced to a combinaison of A, 

B, C, V 


15 




figure 2 


Figure 2: A graphic representation of the relation between the amplitudes A, 
B, C, V when the held is charged. The relation between the amplitudes 
is more complicated than in hgure 1 because of the non vanishing amplitude 
(5^ to produce pairs of /i, p, quanta. The conventions are the same as in 
hgure 1, except that n and p denote a particle and antiparticle of the 
held, and are represented by slanted lines since they are also accelerated by 
the electric held. Note that this picture is not the unique representation of 
the amplitudes A, B, C, V since as indicated in the text each amplitude 
corresponds to 4 diherent processes. We have used this to represent the 
amplitude A{—k\ — k', —k") as the amplitude fh{k') + p{k") M{k), but it 

could equally well be represented by M{—k) m{—k') + jj,{—k''). 
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still vanishes simply because due to charge conservation the quadratic phases 
in t cancel each other and the remaining t integral yields to a delta function 
whose argument is strictly positive on the domain of integration. This iden¬ 
tity implies relations between the amplitudes A, B, C, V which generalize eq. 
Indeed from eqs {24) and (|63l) we obtain immediately: 


V(| - k, k\ k") = l3MHk\k\ k"). (64) 

On the other hand, starting from eq. (^) and using the identity (^4]) to split 
into sums, successively, the functions [A'],P* [A"] and A'] and 

hnally using the complex conjugate version of relation (^), we obtain: 

V(| - k, k', k") = [imB{k'\k,-k") + l3^C{k"\k, -k') - p^P^Ai-k\ - k',-k") (65) 


These relations are illustrated in figure 2. Using eq. (|64D to eliminate V(| 
k, k\ k”) yields 


B(P|fc, -k") _ f3M C(fc"|fc, -k') A(-fc| - P, -k”) 
K{k\k',k'') (3m (3m A{k\kkk") ^ A{k\k'A”) 

This shows that when the charge of the exchanged particle tends to zero, i.e. 
(3fj_ = ratio of the amplitudes A{k\k', k") and B(A;'|fc, —k") 

is still directly related to the amplitudes for the Schwinger process: 


B{k'\k.-k") _ Pm 
A{ k\k'A") (3m 


(67) 


This estimate is proven in Appendix 2 where the a ^ 0 limit is studied in 
detail. It is shown that while each amplitude involving a charged /r-particle 
involves polynomial corrections in a/fi, the ratio of the population of m and 
M species is in thermodynamical equilibrium up only to non-perturbative 
corrections, as expressed in eq. (|67|) . 

Thus the equilibrium distribution of the detector states is given by 


TTiW-/V-C/-r7 


( 68 ) 


Therefore this equilibrium distribution is governed by the hnite 
the horizon area, c.f. the Introduction. 

Upon taking the double limit 5M/M —>• 0,a/Q —>• 0 with 
acceleration 


a = 


{^)E 

(M^) 


change of 
the mean 

( 69 ) 
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fixed, one gets the linearized expression 


P, 


M equil 


P 


= exp 


m equil 


Stt , . E , ^ / 9\ 

- — (Am — a—) + 0{a^) 
a 2a 


+ 0(0,) 


(70) 


governed by a temperature a/27r and an electric potential (= E/2a). 

Thus when the charge of the emitted particle is small enough with re¬ 
spect to its mass square (so as to neglect its vacuum instability), its effect, at 
the linearized level, is to modify Unruh equilibrium by the addition of an a 
electric potential thereby enlarging the relation to thermodynamics in a non 
trivial way 
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5 Appendix 1: Amplitudes of decay due to 
the exchange of charged particle 

For completeness we give here closed forms for the amplitudes A, B, C, V. 
We first discuss the A(/c|fc', k”) amplitude. We have to evaluate the integral 
which in terms of Whittaker’s functions, reads as : 


— J dt 

D-ie^-l/2 


-e {t + 


k' 


k 


MA 


)V2MA 


(71) 




D 


-'je^-l/2 


g-3i7r/4(^ ^ _) /^ 

pa * 


The evaluation of this integral is similar to the pattern followed to obtain 


eq. (|40D . First we split the function 
others functions thanks to the relation 


into two 
and reexpress all the parabolic 


cylinder functions involved in terms of their integral representation 
yields: 


This 


Xi = 


g 4 g 4 g 4 

F(iem + l/2)F(ie^ -|- 1/2) a/2^ 
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Jo 


Z^M^-2H^^-i«V2UIJ(t+^)+m2/2 _L ^-^eM JuVWEit+-^)+iuy2 


e 4 “ e 8 e 


+ e 4 8 e 


QB-' 


(72) 


with k = k' + k". Charge conservation Q = q + a eliminates the qnadratic 
terms in t and the t integration leads to two delta functions : 5{y/^v + 
V aEw =F ^/QEu). The positivity of the u, v and w variables makes only the 
hrst one contributing to the amplitude. After some elementary algebra, we 
obtain for Xi the expression; 


Xi = 


r — p-- 


where: 
X 2 = 


with 

and 

and 


T{iera + l/2)r(ie^ + 1/2) ^J2QE 


\/^ itE- 

g 2 QB qE aE 


(73) 


dvdw e ^ VQBqBcB ^ ( 74 ) 


/o 


\/ QE \fQE 

V Q.E V qE -y qE 2 


= 


= 


qEk-QE k' 

^/QEqEaE 

^ T *7^ 6 jv^ 


D_i^_i +v^e*4f2 -I- _^£: -I-_) 

7 ^ /. 1 . 1 T 

2Fi\1€m + 2 ’ 2 ’ ^ ~ ' 


(75) 

(76) 






V^e *4^2 i?(iem +-, —+-) 

1 . 1 . . qE 

2^1 (^Cm + + -, 1 + ~ f) 

Z Z cxl^ 


(77) 


The way to obtain this last result and others similar is postponed to the 
end of this appendix. Collecting all these results, the amplitude A{k\k',k'') 
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reads: 


K{k\k\k'') = 


_ On/ t--* k' k" •^ \ (-^2 

/— - ITT fc _L£ ^ TTC 1. ( ^ _ 

y^TTC 4 g 2 8 ^2^QE qE otE 2 

(2(3£)i/<(25B)3/-i(2aB)i/->r(!e„ + l/2)r(K„ + 1/2) 


:-)4'‘‘(-)’'"(-)’‘"r(!f+ i)l3 

q q a 2 


( 78 ) 


As expected this result is symmetric with respect to the exchange (m, a, k') 

(p, a, A:"). 

The computation of the B(A:'|A:, —fc") is somewhat different. Starting from 
the general expression (|58|) , we have to evaluate the integral: 


Ji = dt D_ 


icM-l/S 




(79) 






D-ie^-l/2 


k'‘ 


aE 


with k = k' + k". 

After the same transformations as those performed to evaluate Xi, we obtain: 


Ji = 


g 4 i^^Tn~\~€.fn') 


e 4 




_k\2 


T[iem + 1/2) r(A^ + 1/2) a/2^ 






e 4 ” e 8 e 


+ e 4 A«e 8 e 


QB'' 


(80) 


As previously charge conservation implies that the quadratic terms in t can¬ 
cel and the f-integration yields a delta functions 6{VaEw — \f^v ± V QEu). 
However, this time, both terms will contribute to Ji and we are left with the 
following expression: 


Ji = 


-f (cm+e^i) gX 


e 4 


r(iem + l/2)r(ie^ + 1/2) 


TT ^ 3i7r ^ 

g4"Mg 8 ^2 + 


a/^( 


_ k '^ (fc ")2 1 

7^g2\QB oB <^B / 


^/2qE 


y/2aE 


Stt ^ _ 

e 4 A^e 8 jTa 


(81) 


where 


^2 = 


dudwe ” VQBgBcB ^vAb ^vAb ^ 


, y/oE nj QE 


pem-l/2 -ieM-l/2,,,ie;3-l/2 


(82) 
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and: 


J 3 = / dudve ^qe.e^e ^ > 

Jo 

,VQE^, , VqE .■ 


zU 




V aE ^/oE 


u 


(83) 


Note that the last integral JTs can be obtained from the first one J '2 by 
exchanging (m, a) with (/i, a) and {k, k', k”) with (—/c, —k", —k'). In par¬ 
ticular, note the invariance qEk — QEk' 1 —>• —aEk + QEk" = qEk — QEk'. 
We display here the result for 


J 2 = + 


-D 


—iS—^ 




(84) 


^ -I- -, —i^M + 2 ) 2-^1 (*'S -I- -,—ieM + 2 ’ ^ ~ 

Therefore, the B(/c'|A;, —/c") amplitude reads: 


with 
J 4 = 


B(A:'|fc,-fc") = 


('^'1 5(6,1—iciu) f(em—ieivf) 


(85) 


o-'|(em+e^)g'T^27re 2 (qe ^E+qe) 


{2QEY/i{2qEY/^{2aEYl^T{iem + l/2)r(ze^ + 1/2) 


. 1 , Tfg _ 

r(i£: -h -)e 2 e « e « 


v^e* O 


1^4 


r_ stTT ,1 i,, ii ^.E , 

|e 8 + -j-icM + -) 2 ^i(^<^ 2’ 2 ~ + 

_1_ g ■"■eMg g B(i€m + 2> + 2 ) 2^1 (*'^ 2 ’ 2 ~ ^ ~ q^)'|(86) 


Finally, notice that the C{k"\k, —k') amplitude can be easily obtained from 
B(/c'|A;, —k") by the substitution m, q, k' h-> /i, a, k'\ and hence hi i—> —hi. To 
complete the calculation we have to evaluate the integrals X 2 , J 2 and As 
an example we consider the integral X 2 , the others follow a similar pattern : 


X 2 = / dvdw e ^ ^QEqEaE ^,/We 

Jo 




vw 

By using the reduced variables: 


l/2„,26,71 1/2 


X = 


/ aX^ 




kqE-QEk' 


v,Y = ^ —w,n= ) ^ , c = = 

’ \/ iO 7-1 ’ /iO 7-1 . 7-1 _ f=i 5 




YQEqEoE 


I qE 

aE 


(87) 


( 88 ) 
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this integral becomes: 


qE 


2e^ + l/2 


IQE 

aE 


iCr7i + l/2 


I 


( 89 ) 


where: 


1= r dXdY + c-iF)-*^"-V2x*^--V2yi^^-i/2_('go) 

Jo 

The new changes of variables: 

Q _ Q 

X = Y = ^ ^ 0 < p < oo, —oo < 9 < oo (91) 


and z = psinh^ factorizes the double integral: 

1 '91 1 

J = 2-*^+2 / dzd 9 e- 2 ^ z^^-^{smh 9)-^^-2 

Jo 

|g*6»(e^-e;,)j'g0+7 _|_ ^-O+'y-s^-ieM-^ 

^ g-0+7^-ieM-5g*\/2O2;| 


( 92 ) 

( 93 ) 


where the two terms come from the separation of the positive and negative 
values of the integration variable 6. The z integrals are the representation 
(1^) of the parabolic cylinder function, so: 


O T 2 /* OO 

I = 2r(i-+/ d0(2sinh0)-*^-^ (94) 

2 2 Jo 


A hnal change of variable s = reduces the last integrals to products of 
Euler and hypergeometric functions 
for I 2 reads: 


211, formula 3.197.2 ). The final result 


I 2 = 


I QE 

qE 


ie^+ 1/2 


QE 

aE 


t + 1/2 


if 1 

riy + 2)6 "6 4 6 - 


(95) 


+ D-i£-L 


-V2E^n 


111 . 1 . aE 

B(ie^ + —5 "I" o) 2-^1 (7: + '^^fi + —, 1 — icm + i^M] -77) 

Z Z Z Z qh/ 


/. 1 lx . . 1 . 

Byiem + 7:5 + 7:) 2-^i(7: + + 7:, 1 — - 

Z Z Z Z ah 


Similarly, for the amplitude B(fc'|fc, —k") the last integration can be carried 
out by using (0, formula 3.197.1). 
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6 Appendix 2: Small charge limit of the am¬ 
plitudes 

In this appendix we compute the small charge limit of the amplitudes A{k\k\ k") 
and — fc"). We also prove eq. (|67D , i.e. that the ratio of the ampli¬ 
tudes A{k\k\k") and —fc") is given by the ratios of the Schwinger 

amplitudes, up to terms proportional to Moreover we shall evaluate 

A{k\k'^k") in the limit a —>■ 0. As emphasize in the main text, a consequence 
of this calculation is that while the amplitudes A(/c|fc', k") and B(/c'|fc, —k") 
differ polynomially in the variable a/fi from their chargeless limits, their ra¬ 
tio is nevertheless given by eq.(^), which differs from the chargeless limit 
only by non-perturbative corrections. When a is small, the two hypergeo¬ 
metric functions giving J '4 in eq (^) can be combined into one thanks to the 
formulae (9.131.2 and 9.131.1)of ref. 1^ : 


gl7r/8g W€M y/ctE , -— 

J> = --(n/^) 




QE)' 


•B(, 




-2Fi{- + iem,- + teMA 


+ e~ 




+ Fm, 


Note that (and this constitutes a check of the exactness of the calculation), 
that in the limit a/fi ^ 0 , there appears diverging phases in the hrst factor 
of this expression which cancel each other. If we omit the small corrections 
proportional to the Schwinger factor most of the prefactors are com¬ 

mon between the amplitudes A{k\k'k'') and B{k'\k — k") and thus can be 
omitted in the ratio A/B which reduces to: 

A{k\k'k") ^ T{\-i£)T{\ + i£)e^M 


B{k'\k-k") 


+ e" 


'OA/fi) 


V{ie^ - zeM)r(l - ie^ + ieM) 

= + e-^^^0{a/fi) (97) 

which is the sough for result. Now we discuss in more details the limit 
of the amplitude A{k\k',k") when a —>■ 0, and check that we recover the 
a = 0 result in the limit. The computation is done in three steps. Eqs 
(|78|j77|) express A{k\k',k'') as products of phases and T functions with a 
sum of products of Eulerian (B) functions, Whittaker’s (V) functions and 
hypergeometric functions. 

Firstly we obtain, by a saddle point evaluation, appropriate approximations 
of the Whittaker’s functions. Then we estimate the limits of the hyperge¬ 
ometric functions as confluent hypergeometric functions. Finally, applying 


qE 

aE’ 

(96) 
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several times the Stirling and reflection formnlae on the F- and 5-Enler 
fnnction we obtain the resnlt. 

The integrals 


r(y + ^)e^e^e-^^D_.£_i [±V2e^^n\ = dv 


have saddle points located respectively at: 


+ £)l/2 

Vs = - 

2 

They are approximated by 







(99) 


( 100 ) 


with 

I 3 (fc + fc')^ (fc + fc')w £ w =F (fc + fc') (k + k')(k — k') (k — k')oj . . 

4 aE AaE 2aE 2 ^2^ 2QE 2QE ^ ’ 


Let ns emphasize that both integrals (p8D have their modnlns of the same 
order of magnitnde in a//i. 

To obtain the limits of the hypergeometric fnnctions appearing in eq.(^7D 
as continent hypergeometric fnnctions when a —>■ 0 is straightforward. For 
small valne of a we obtain : 


1 1 ctE 1 Q; 

2Fi(ieM + 77) i^fj. + 77) 1 + -77) = M(icM + 77) 1 + — *em) —77—77) + O(~)(102) 

2 2 qE 2 2qE fi 


and: 


„ / . 1 . 1 1 . qE . -i-nLr + , 

2J^i(jeM+2)*em + -)l + ieM-«eM;-^) = e ^(^^) " 


(eM+Crr 


w 


pem+eM),'J(em-£M) 


(103) 


+ o(^^). 


Here also these two fnnctions are of the same order of magnitnde in a/ju 
bnt the prefactors mnltiplying them in the amplitnde A(kjk', k") are qnite 
different. The first one is mnltiplied by: 


+ 0 > + *(«M - em) + 77 ) = 


27j-gi(eM-em)g7reM/2g-7rem/2g-7re^ 

r(l + i(eM ~ £m)) 


l + 0(-) (104) 
A* 


whereas the second one is pondered by: 

B(iern + i + ^(eM - Cm) + i) = r(*e™ + i)e--^-i/2g—„/2gW4(i + O(“)).(l05) 

Z A Z fJj 
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We see that the exponential factor appearing in eq.( [104|) makes the hrst 
term of X 3 in eq.(^) negligible with respect to the second one in the limit 
a —> 0 at hxed non vanishing valne of /r. Collecting all the results (| 102 | , |1U4| ) 
and (1^) , we obtain the limit we are discussing. At zero order in a/zr, once 
more (as expected) all the diverging phases cancel in the hrst term and the 
remaining factors group together to give the expression (^) with qE = QE. 
So, we obtain at the end: 


A{k\k'k'') = {Ao{k'\kk'') + 0{a/fi)) + (106) 

where Ao{k\k'k") is the transition amplitude in the neutral case (^Ol). Using 
eq. (|56|), one can similarly show that 

B{k\k'k'') = {Bo{k'\kk'') + 0(a/zi)) + (107) 


where Bo{k\k'k'') is the transition amplitude in the neutral case, from these 
relations (|106|, p.07|) one can only a priori deduces that their ratio behaves as: 


A{k\k'k") 

B{E\k-E') 


g—+C)(-)) + e-"^'' 0 (a/zi) . (108) 

/i 


while, our previous computation, eq. (|^ , shows that actually all the poly¬ 
nomial corrections in a/p to the boltzmannian factor cancel each others! 
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